Singlet Excitations in Pyrochlore: A Study of Quantum Frustration 
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We apply the Contractor Renormalization (CORE) method to the spin half Heisenberg an- 
tiferromagnet on the frustrated checkerboard and pyrochlore lattices. Their ground states are 
spin-gapped singlets which break lattice symmetry. Their effective Hamiltonians describe fluc- 
tuations of orthogonal singlet pairs on tetrahedral blocks, at an emergent low energy scale. We 
discuss low temperature thermodynamics and new interpretations of finite size numerical data. 
We argue that our results are common to many models of quantum frustration. 
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Frustration in classical spin models often leads to a 
complex energy landscape. Certain models, such as the 
Heisenberg antiferromagnet on the pyrochlore lattice, has 
an extensively degenerate ground state manifold. This 
model given by 



(i) 



has spins sitting on corner sharing tetrahedral units 
(see Fig. 0(a)). 

In the semiclassical approximation |2| [3|, a large de- 
generacy survives the quantum fluctuations, and thus re- 
sists ground state selection by the "order from disorder" 
mechanism. 

A pressing open question is what happens in the 
strong quantum limit, e.g. the spin half case? Series 
expansions^ suggest rapid decay of spin correlations. 
Does this indicate the formation of a translationally in- 
variant spin liquid or lattice symmetry breaking valence 
bond solid? In the quantum case, is there an emerging 
low energy scale, in lieu of the classical ground state de- 
generacy? 

The purpose of this Letter is to derive the low en- 
ergy effective Hamiltonian starting from the Heisenberg 
model. As a warm-up to the pyrochlore lattice (Fig. 
0(a)) , we treat its two dimensional (toy model) re- 
duction, the Checkerboard lattice (Fig. 0(b)). The 
Checkerboard has recently received significant theoret- 
ical attention||, 0, g, g 

Our approach is the Contractor Renormalization 
(CORE) methodp). The CORE is a real-space discrete 
renormalization transformation invented by Morningstar 
and Weinstein. It maps a lattice Hamiltonian to an ef- 
fective Hamiltonian with the same low energy spectrum. 
The CORE method computes the effective interactions 
at all ranges using exact diagonalizations of finite con- 
nected clusters. Truncation of interactions beyond a fi- 
nite range is an approximation whose error can be es- 
timated numerically from the next higher range terms. 
CORE has been successfully applied to describe the spec- 
tra of Heisenberg models on chains and ladders JO], [l2| . 



Recently, it was applied to the square lattice Hubbard 
model to derive the Pla quet te Boson-Fermion Model for 
cuprate superconductors Jl3| . We refer the reader to pre- 
vious reviews jll], [l3| for the mathematical background 
and technical details. 

For each of the Hamiltonians at hand, we define local 
operators from the lowest eigenstates of the elementary 
clusters, e.g. a tetrahedral unit in the Pyrochlore lattice. 
We shall compute the effective interactions by CORE up 
to four clusters range, and estimate the truncation error. 




(a) (b) 

FIG. 1: The Pyrochlore (a) and Checkerboard (b) lattices. 

Our key results are the following: For the Checker- 
board, we confirm the conclusions of refs 
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by other means, that the ground state is closely approx- 
imated by the product of uncrossed plaquette singlets. 
The effective pseudospin Hamiltonian allows us to inter- 
pret the numerical spectrum of low lying singlets^], 0] 
in terms of Ising domain walls. This yields the number 
of singlets as a function of lattice size: a numerically 
testable prediction. In addition we expect a branch of 
weakly dispersive triplets at an energy scale slightly be- 
low the Heisenberg exchange. For the Pyrochlore, we 
apply two CORE steps to arrive at an Ising-like model 
of local singlets on larger super-tetrahedra blocks which 
form a cubic superlattice. At mean field level, we obtain 
a singlet ground state which breaks lattice symmetry as 
depicted in Fig.||. Here too, the effective Hamiltonian 
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describes Ising-like domain walls. We shall conclude that 
lattice symmetry breaking and local singlet excitations 
are general features of highly frustrated quantum anti- 
ferromagnets. 

The Checkerboard. The lattice depicted in Fig. |l|(b), 
contains crossed plaquettes (which are equivalent to three 
dimensional tetrahedra) connected by uncrossed plaque- 
ttes. The first step is to choose elementary clusters which 
cover the lattice. Two different options for plaquettes are 
the crossed and the uncrossed plaquettes. While the two 
clustering choices may appear to yield different ground 
states and excitations, we shall see that they are in fact 
consistent with each other, and yield complementary in- 
formation. 

Clustering with uncrossed plaquettes. From the spec- 
trum of a single uncrossed plaquette, we retain the sin- 
glet ground state as a vacuum state | } i and the low- 
est triplet as a singly occupied boson state t' ai | O ) i . 
a — x, y, z is a cartesian index of the triplet. 

The effective Hamiltonian in the uncrossed plaquettes 
basis is (in units of J = 1): 

(ij)a/3 

(2) 

where Lj = Ylaff^aj^otptpj and £ap are 3 x 3 spin-1 
matrices in a cartesian basis. 

The parameters calculated by CORE upto range 2 are: 
e' t = 0.5940, K = 0.2985, A = 0.1656 and b = 0.0776. 
The truncation errors from up to range-4 are less than 
2%, and will be ignored Q. 

Note that H-eff in (@) commutes with the number of 
triplets since it has no anomalous pair creation terms, as 
appear e.g. for the square lattice JD|. Thus, at this level 
of truncation, the plaquette vacua product 

i*o>=ni n <>> ( 3 ) 

i 

is an exact ground state of the effective Hamiltonian 
(|)(T|. This result agrees with Moessner et al ||, who 
argued for a plaquettized singlet ground state based on 
an effective quantum dimer model. 

We are also able to obtain the triplet (spin) gap for 
Since i |0) is an approximate eigenstate of (|J), its en- 
ergy (spin gap) can be read from e' t — 0.5940. This com- 
pares well with the value of 0.6-0.7 estimated by exact 
diagonalizations of finite systems fit] . We have found very 
weak hopping terms (of magnitude 0.01J) due to CORE 
interactions of range four, which will give the triplets a 
weak dispersion in the full lattice. 

Clustering with crossed plaquettes. The isolated 
crossed plaquette has two fold degenerate singlet ground 
states, which we can represent by a pseudospin-i doublet 
(see Fig. |). 



The quantization axis for the pseudospin operators is 
chosen as in Ref. Jll|, with the +z(—z) directions rep- 
resenting states with positive (negative) chirality. The 
planar angles 0, 7r/3, 27r/3, represent the three (non or- 
thogonal) dimer configurations of the tetrahedron. The 
states with their pseudospin polarized in the +x and —x 
direction are shown in Fig. |2[ 

The effective Hamiltonian in the crossed plaquettes ba- 
sis is an Ising-like model 

Heff = -J' ■ «Hi)(Sj ■ e v) - h E S i- ( 4 ) 

(y) « 

where are directors on the x — y plane pointing 
7r/3 (— 7t/3) away from the x axis for horizontal (verti- 
cal) bonds. At range-2 we obtain, J' = 0.527665, and 
,h = 0.118084 Corrections from range-3 and range-4 
CORE were computed Q|, and found to be unimportant 
for the symmetry and correlations of the low excitations. 

We can solve ([|) in mean field theory. The energy ex- 
hibits two minima, where the pseudospins describe ver- 
tical or horizontal dimers. These states correspond re- 
spectively to projections of the two equivalent plaquette 
ground states onto the truncated Hilbert space of the 
crossed plaquettes. Although the ground state energy is 
not well converged at range-2 CORE, (||) treats the two 
symmetry breaking ground states in an unbiased fash- 
ion. It therefore describes low energy singlet excitations 
which are pseudospin-flips or Ising domain walls between 
ground states (see Fig. [|). 

Thus we see that the two methods elucidate compli- 
mentary aspects of the Checkerboard. The approach us- 
ing uncrossed plaquettes gives a very accurate descrip- 
tion of one ground state, which serves as its vacuum, but 
where the other ground state is a multi-magnon bound 
state. In contrast, the crossed plaquettes describes cor- 
rectly the low energy singlet excitations by an effec- 
tive Ising model J Ismg S z S z with a coupling constant 
jismg _ j' renormalized by the quantum fluctua- 
tions. The lowest excitations are gapped spin flips of 
energy J Ism 9 , They gain a weak dispersion due to the 
effective S + S~ couplings. 
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FIG. 2: The two singlet ground states of the uncrossed pla- 
quette and their pseudospin representation. In the physical 
state thick lines denote valence bond singlets. 
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FIG. 3: A singlet excitation of the checkerboard. This state 
is a domain wall (dashed line) between domains of the two 
plaquettized ground states. The corresponding pseudospin 
directions in the crossed plaquette approach are depicted by 
arrows. 

The number of singlet states is expected to grow as 
power laws with the size of the system N. For exam- 
ple, the lowest lying single spin-flips grow as N, while 
higher spin flip pairs grow as TV 2 and so on. An Ising- 
like phase transition is expected between the broken and 
unbroken translational symmetry phases at a tempera- 
ture scale of T c ~ J Ism 9 1 with an associated logarithmic 
divergence of the heat capacity at T c . 

The Pyrochlore, depicted in Fig. [l](a), is a three di- 
mensional network of corner sharing tetrahedra. Like the 
checkerboard, it has a macroscopically degenerate (expo- 
nential in lattice size) classical ground state manifold. 
For the quantum S=l/2 case, local pseudospins can be 
defined by the degenerate singlets of disjoint tetrahedra. 
These cover all sites of the Pyrochlore and form an FCC 
superlattice. The effective hamiltonian on this FCC lat- 
tice was calculated by CORE. The first non-trivial inter- 
teterahedra coupling are obtained at range three con- 
nected tetrahedra, which yield 

H e% =£ ((^(S,-egi)(S i .eg)+ (5) 

(ijk) \ 




The coupling parameters (in units of J) are: J2 = 
0.1049, J 3 = 0.4215, and e$ 3 ,i = 1,2,3 are three unit 

(i) 

vectors in the x-y plane whose angles cv^ 23 depend on 
the particular plane defined by the triangle of tetrahc- 
dral units 123 as given in table I of |l6| . The effective 
hamiltonian (|5|) resembles the terms obtained by second 
order perturbation theory (in inter-tetrahedra couplings) 
by Harris et al [Q and Tsunctugu[l6|. The classical 
mean field ground state of (||) is identical to the ground 
state found in Ref . p6{ : three of the four FCC sublattices 
are ordered in the directions e(0), e(27r/3), e(— 27r/3), 
while the direction of the fourth is completely degener- 
ate. Therefore, classical mean field approximation for 



(g) is insufficient to remove the ground state degeneracy. 
Tsunetsugu|ll| was able to lift the degeneracy by includ- 
ing spinwave fluctuations effects which produce ordering 
at a new low energy scale. 

Here we avoid the a-priori symmetry breaking needed 
for semiclassical spinwave theory, by treating ^) fully 
quantum mechanically. This entails a second CORE 
transformation which involves choosing the "supertetra- 
hedron", as a basic cluster of four tetrahedra, whose struc- 
ture and spectrum are depicted in Fig.||. 

Our new pseudospins Tj are defined by the two degere- 
nate singlet ground states of the supertetrahedron. (This 
degeneracy is found for the Hciscnbcrg model on the orig- 
inal lattice as well as for the effective model (||)). These 
states transform as the E irreducible representation of 
the tetrahedron (Td) symmetry group, similarly to the 
singlet ground states of a single tetrahedron. 
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FIG. 4: Low energy spectrum of a single "supertetrahedron". 
A dashed line illustrates the energy cutoff: only the two de- 
generate ground states are retained in the Hilbert space. 

The supertetrahedra form a cubic lattice, shown in Fig. 
||. The effective hamiltonian (||) and the lattice geome- 
try imply that non-trivial effective interactions appear 
only at the range of three supertetrahedra and higher. 
Range three effective interactions include two and three 
pseudospin interactions, which are dominated by 

We// = Ji£(n-%)(r J --f ij )+ (6) 

(ij) 
{(ij)) 

J 2 ] £ X ' fly x ' x *))• 

Here, ( ) and (( }} indicate summation over nearest- and 
next nearest-neighbors, respectively. The coupling con- 
stants are found to be relatively small: J\ = 0.048J, 
j( a) = -0.006J and = 0.018J . The vectors fy de- 
pend on the vector connecting the two sites, and their 
values are presented in table . 

We performed classical Monte Carlo simulations us- 
ing the classical (large spin) approximation to (||). The 
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TABLE I: The values of the vectors fij in eq. (|6|), depending 
on the vector r^- separating the sites i and j. 



(±1,0,0) 
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FIG. 5: The pyrochlore viewed as a cubic lattice of supertetra- 
hedra. The arrows show the direction of the supertetrahedra's 
pseudospin in the MF ground state. 



ground state was found to choose an antiferromagnetic 
axis, and to be ferromagnetic in the planes as depicted 
in Fig. ||. This ground state differs from the semiclas- 
sical ground state The latter involves condensa- 

tion of high energy states of the supertetrahedron in the 
thermodynamic ground state. Since in the supertetrahe- 
dra diagonalizations, we find a much larger energy gap 
than inter-site coupling we believe these excitations can- 
not condense to yield the semiclassical ground state sym- 
metry breaking. 

Discussion. The CORE technique enabled us to derive 
effective Hamiltonian for highly frustrated antiferromag- 
nets, written in terms of low energy, local degrees of free- 
dom. For both Checkerboard and Pyrochlore systems, 
we found lattice symmetry breaking ground states which 
are essentially products of local singlets. The spin gap 
to the lowest triplet excitation is large (of order J), and 
seems to survive interplaquette interactions. The low en- 
ergy excitations are singlets, which are local pseudospin 
flips, or Ising domain walls between ground states. The 
ordering energy scale is of order J/100 for the Pyrochlore. 

This picture seems to be consistent with existing nu- 
merical data for the density of low energy singlets on the 
Checkerboard and Pyrochlore. Experimentally, lattice 



symmetry breaking could drive a static lattice distortion, 
which would be observable by additional Bragg peaks in 
neutron and X-ray scattering. For example, the anti- 
ferromagnetic order between planes of supertetrahcdra 
would correspond to a lattice distortion with wavelength 
of four tetrahedra. 

How general are these results ? Formation of local sin- 
glets is a natural way to relieve the frustration in quan- 
tum antiferromagnets that can be written as a sum over 
clusters 2 c (£i Sic) 2 - On each even cluster the ground 
state is a singlet with a gap of order 1/5 to a local 
triplet. Frustration suppresses hopping of these triplets 
and could inhibit their condensation into a spin ordered 
ground state. Thus lattice symmetry breaking singlet 
ground states are expected as a typical feature of frus- 
trated quantum antiferromagnets |17[ . 
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